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Abstract 

In the supersymmetric extension of the Standard Model we derive the two-loop 
QCD corrections to the scalar quark contributions to the electroweak precision ob- 
servables entering via the p parameter. A very compact expression is derived for the 
gluon-exchange contribution. The complete analytic result for the gluino-exchange 
contribution is very lengthy; we give expressions for several limiting cases that were 
derived from the general result. The two-loop corrections, generally of the order of 10 
to 30% of the one-loop contributions, can be very significant. Contrary to the Standard 
Model case, where the QCD corrections are negative and screen the one-loop value, the 
corresponding corrections in the supersymmetric case are in general positive, therefore 
increasing the sensitivity in the search for scalar quarks through their virtual effects in 
high-precision electroweak observables. 
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1. Introduction 



Supersymmetric theories (SUSY) |jl|] are widely considered as the theoretically most ap- 
pealing extension of the Standard Model (SM). They are consistent with the approximate 
unification of the three gauge coupling constants at the GUT scale and provide a way to 
cancel the quadratic divergences in the Higgs sector hence stabilizing the huge hierarchy 
between the GUT and the Fermi scales. Furthermore, in SUSY theories the breaking of the 
electroweak symmetry is naturally induced at the Fermi scale, and the lightest supersym- 
metric particle can be neutral, weakly interacting and absolutely stable, providing therefore 
a natural solution for the Dark Matter problem; for recent reviews see for instance Ref. 0. 

Supersymmetry predicts the existence of scalar partners fi, fn to each SM chiral fermion, 
and spin-1/2 partners to the gauge bosons and to the scalar Higgs bosons. So far, the direct 
search of SUSY particles at present colliders has not been successful. One can only set lower 
bounds of C(IOO) GeV on their masses 0. The search for SUSY particles can be extended 
to slightly larger values in the next runs at LEP2 and at the upgraded Tevatron. To sweep 
the entire mass range for the SUSY particles, which from naturalness arguments is expected 
not to be larger than the TeV scale, the higher energy hadron or e^e~ colliders of the next 
decade will be required. 

An alternative way to probe SUSY is to search for the virtual effects of the additional 
particles. Indeed, now that the top-quark mass is known g], and its measured value is 
in remarkable agreement with the one indirectly obtained from high-precision electroweak 
data, one can use the available data to search for the quantum effects of the SUSY parti- 
cles: sleptons, squarks, gluinos and charginos/neutralinos. In the Minimal Supersymmetric 
Standard Model (MSSM) there are three main possibilities for the virtual effects of SUSY 
particles to be large enough to be detected in present experiments: 

i) In the rare decay b 57, besides the SM top/W-boson loop contribution, one has 
additional contributions from chargino/stop and charged Higgs/stop loops 0. These 
contributions can be sizable but the two new contributions can interfere destructively 
in large areas of the MSSM parameter space, leading in this case to a small correction 
to the decay rate predicted by the SM. 

ii) If charginos and scalar top quarks are light enough, they can affect the partial decay 
width of the Z boson into b quarks in a sizable way [Q. This feature has been widely 
discussed in the recent years, in view of the deviation of the Z —>■ bb partial width 
from the SM prediction 0. However, for chargino and stop masses beyond the LEP2 
or Tevatron reach, these effects become too small to be observable 

iii) A third possibility is the contribution of the scalar quark loops, in particular stop and 
sbottom loops, to the electroweak gauge-boson self-energies [H, |^: if there is a large 
splitting between the masses of these particles, the contribution will grow with the 
square of the mass of the heaviest scalar quark and can be very large. This is similar 
to the SM case where the top/bottom weak isodoublet generates a quantum correction 
that grows as the top-quark mass squared. 
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In this paper, we will focus on the third possibility and discuss in detail the leading con- 
tribution of scalar quark loops to electroweak precision observables, which is parameterized 
by their contribution to the p parameter. The radiative corrections affecting the vector boson 
self-energies stemming from charginos, neutralinos and Higgs bosons have been discussed in 
several papers @, 0. In the MSSM, because of the strong constraints on the Higgs sector, 
the propagator corrections due to Higgs particles are very close to those of the SM for a light 
Higgs boson []TU|. In the decoupling regime where all scalar Higgs bosons but the lightest 
are very heavy, the SUSY Higgs sector is effectively equivalent to the SM Higgs sector with 
a Higgs-boson mass of the order of 100 GeV. The contribution of charginos and neutralinos, 
except from threshold effects, is also very small p. The main reason is that the custodial 
symmetry which guarantees that p = 1 at the tree level is only weakly broken in this sector 
since the terms which can break this symmetry in the chargino/neutralino mass matrices 
are all proportional to Mw and hence bounded in magnitude §]. 

The propagator corrections from squark loops to the electroweak observables can be 



attributed, to a large extent, to the correction to the p parameter [|lT|], which measures the 
relative strength of the neutral to charged current processes at zero momentum-transfer. 
This is similar to the SM, where the top/bottom contribution to the precision observables 
is, to a very good approximation, proportional to their contribution to the deviation of the p 
parameter from unity. Further contributions, compared to the previous one, are suppressed 
by powers of the heavy masses. It is mainly from this contribution that the top-quark mass 
has been successfully predicted from the measurement of the Z-boson observables and of the 
W-boson mass at hadron colliders. However, in order for the predicted value to agree with 



the experimental one, higher-order radiative corrections [0, 13, had to be included. 
For instance, the two-loop QCD corrections lead to a decrease of the one-loop result by 
approximately 10% and shift the top-quark mass upwards by an amount of ~ 10 GeV. 

In order to treat the SUSY loop contributions to the electroweak observables at the 
same level of accuracy as the standard contribution, higher-order corrections should be 
incorporated. In particular the QCD corrections, which because of the large value of the 
strong coupling constant can be rather important, must be known. In a short letter |]l5l we 
have recently presented the results for the 0{as) correction to the contribution of the scalar 
top and bottom quark loops to the p parameter. In this article we give the main details of 
the calculation and present the explicit result for the gluon-exchange contribution as well 
as the result for the gluino-exchange contribution in several limiting cases. 

The paper is organized as follows. In the next section, we summarize the one-loop results 
and fix the notation. The main features of the two-loop calculation are discussed in section 3. 
In section 4 a compact expression is given for the gluon-exchange contributions. The results 
for the gluino-exchange contributions are presented for the limiting cases of zero gluino 
mass and a very heavy gluino as well as for the case of arbitrary gluino mass but vanishing 
squark mixing. Effects of 0{as) corrections to relations between the squark masses existing 
in different scenarios are discussed. In section 5 we give our conclusions. 
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2. One- loop results 

For the sake of completeness, we summarize in this section the one-loop contribution of a 
squark doublet to the electroweak precision observables. Before that, to set the notation, we 
first discuss the masses and couplings of scalar quarks in the MSSM. 

As mentioned previously, SUSY associates a left- and a right-handed scalar partner to 
each SM quark. The current eigenstates, ql and (Jr, mix to give the mass eigenstates qi and 
q2', the mixing angle is proportional to the quark mass and is therefore important only in the 
case of the third generation squarks. In the MSSM, the squark masses are given in terms 
of the Higgs-higgsino mass parameter /i, the ratio of the vacuum expectation values tan/3 
of the two Higgs doublet MSSM fields needed to break the electroweak symmetry, the left- 
and right-handed scalar masses Mq^ and Mq^ and the soft-SUSY breaking trilinear coupling 
Aq. The top and bottom squark mass eigenstates and their mixing angles are determined 
by diagonalizing the following mass matrices 

■'^'^ ~ \ mt Mt"" Ml + m? + I cos 2/3 s^, M| j 



■^^^ ~ 1^ rrif, Mt"" M?^ + < - I cos 2/3 sl^ Ml j ' 

with M^^ = At — fi cot (3 and M^^ = Ai, — fi tan /? ; = 1 — = sin^ 6w Furthermore, 
SU(2) gauge invariance requires M^^ = M^^ at the tree-level. Expressed in terms of the 
squark masses m^^ , irtq^ and the mixing angle 9q the squark mass matrices read 

, _ / cos^ eqml + sin^ ^,m| sin 6, cos ^,(m?^ - m|) \ 
- \ sin^,-cos^,-(m?^ -m?J sin^ ^.m^^ + cos^ ^,-m| ) ' 

Due to the large value of the top-quark mass mt, the mixing between the left- and 
right-handed top squarks II and tj? can be very large, and after diagonalization of the mass 
matrix the lightest scalar top-quark mass eigenstate ti can be much lighter than the top 
quark and all the scalar partners of the light quarks . The mixing in the sbottom sector is 
in general rather small, except if A^y, fi or tan /3 are extremely large. In most of our discussion 
we will assume that, because of the small bottom mass, the mixing in the sbottom sector is 
negligible and therefore Bl = &i- 

Using the notation of the first generation, the contribution of a squark doublet u, d to the 
self-energy of a vector boson V = Z,W and to the Z-j mixing is given by the diagrams 
of Fig. |I|. Summing over all possible flavors and helicities, the squark contribution to the 
transverse parts of the gauge-boson self-energies at arbitrary momentum-transfer can be 
written, in terms of the Fermi constant Gp, as follows: 
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with the reduced couphngs of the squarks to the W and Z bosons, including mixing 9q 
between left- and right-handed squarks, given by (cg and are the electric charge and the 
weak isospin of the partner quark) 

_ I cos 9u cos 9^ — cos 9u sin 9^ \ 
9wudj - I -cos^ssin^j sin sin ) ' 
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(5) 



In both the dimensional regularization [|T^ and dimensional reduction [|T8l schemes,|_] the 
function no(g^, m^,m^) reads 



no(g ,m^,mb) = - 



2 , _2 <f , (2 , 2 <f iK-^ir 



ma + m^, - — + {m^ + m^ - — 



2g2 



ml -ml 



2g2 



{Ao{ma) - Ao(mb)) - Ao{ma) - A^im^) 



(6) 



The Passarino-Veltman one- and two-point functions are defined as 
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Bo{q , m„, mfe) = - + ^^'"(g , m„, mt) + e Bl{q , m^, mb) 
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(3^/^ {q^, ml, ml) ml + ml - q^ + (3^^'^{q^, ml, ml) 
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(7) 



"'^In general the dimensional reduction scheme, which preserves SUSY, should be used. For all quantities 
considered in this paper dimensional reduction yields precisely the same result as dimensional regularization 
(see however the discussion in section 4.4). 
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where /i is the renormahzation scale, [3 the phase space function, 



n/ 2 2 2\ 2 2 2 , ' m'i}'^ 

(3{q , m^)=q -m^-m^,+ 



and 2e = 4 — n with n the space-time dimension. We have absorbed a factor (e'^/47r)'^, 
with 7 the Euler constant, in the 't Hooft scale /i to prevent uninteresting combinations of 
ln47r, 7 ... in our results. The explicit form of the function Bq is not needed for our purposes 
but can be found in Ref. . 

At zero momentum-transfer, = 0, the function Ho is finite and reduces to 

Foiml, ml) = no(0, ml ml) =ml + ml- in (9) 

— mi mi 

We now focus on the contribution of a squark doublet to the p parameter. It is well- 
known that the deviation of the p parameter from unity parameterizes the leading universal 
corrections induced by heavy fields in electroweak amplitudes. It is due to a mass splitting 
between the fields in an isospin doublet. Compared to this correction all additional con- 
tributions are suppressed. In the relevant cases of the W-boson mass and of the effective 
weak mixing angle sin^ 6'^, for example, a doublet of heavy squarks would induce shifts 
proportional to its contribution to p, 

5Mv,^^^4^Ap; 5sin^C--^^Ap. (10) 

In terms of the transverse parts of the W- and Z-boson self-energies at zero momentum- 
transfer, the squark loop contribution to the p parameter is given by 

1 . nzz(o) n,„„(o) ,,,, 
"'—p- ^'^-Si Ml-- 

Using the previous expressions for the W- and Z-boson self-energies and neglecting the 
mixing in the sbottom sector, one obtains for the contribution of the tjh doublet at one-loop 
order (only the left-handed sbottom 6^ contributes for = 0) : 



Apo^usY ^ ^ _sin^^,cos^^,Fo(m|,m|) + cos^^,Fo(m| 



2^ , .^.2a ^,^2 2 



L ' 



+ sin2^,-Fo(m|,m?j] . (12) 



The function Fq vanishes when the two squarks running in the loop are degenerate in mass 
il,ml 



Fo(m^,m^) = 0. In the limit of large squark mass splitting it becomes proportional to the 



heavy squark mass squared: Fo(m^, 0) = m^. Therefore, the contribution of a squark doublet 
becomes in principle very large when the mass splitting between squarks is large. This is 
exactly the same situation as in the case of the SM where the top/bottom contribution to 



the p parameter at one-loop order, keeping both the t and b quark masses, reads |Tl 



Apr = ^Joimlml). (13) 
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For rrit ^ this leads to the well-known quadratic correction Ap^^ = 3G'i7m^/(8v^7r 

In Fig. ^ we display the one-loop correction to the p parameter that is induced by 
the i/b isodoublet. The scalar mass parameters are assumed to be equal, Mj^ = M^^ 
and M^^ = M-^^, as it is approximately the case in Supergravity models with scalar mass 
unification at the GUT scale |]21]]. As mentioned above, SU(2) gauge invariance requires 
at the tree-level[] M^^ = M^^ , yielding in this case M^^ = M^^ = M^^ = M^^ = rriq. In 
this scenario, the scalar top mixing angle is either very small, 9i ~ 0, or almost maximal, 
6i ~ — 7r/4, in most of the MSSM parameter space. The contribution Apg^^"^ is shown as 
a function of the common squark mass iriq for tan/? = 1.6 for the two cases M^^ = (no 
mixing) and M^"^ = 200 GeV (maximal mixing) ;0 the bottom mass and therefore the mixing 
in the sbottom sector are neglected leading to m^^ = m^^ ~ rriq. Here and in all the numerical 
calculations in this paper we use nit = 175 GeV, Mz = 91.187 GeV, M^y = 80.33 GeV, and 
as = 0.12. The electroweak mixing angle is defined from the ratio of the vector boson masses: 

As can be seen, the correction is rather large for small rrig, exceeding the level Ap = 
1.3 X 10~^ of experimental sensitivity^ in the case of no mixing for rriq ~ 150 GeV, which 
corresponds to the experimental lower bound on the common squark mass 0], and getting 
very close to it in the case of maximal mixing. For large rriq values, the two stop and the 
sbottom masses are approximately degenerate since rriq ^ rrit, and the contribution to Ap 
becomes very small. 

For illustration, we have chosen the value tan/3 =1.6 which is favored by 6-r Yukawa 



coupling unification scenarios In fact, the analysis depends only marginally on tan/3 
if M/"^ (and not At and p) is used as input parameter. The only effect of varying tan (3 is 
then to slightly alter the D-terms in the mass matrices eqs. (1-2), which does not change the 
situation in a significant way. However, for large tan/3 values, tan/3 ~ rrit/mb, the mixing 
in the sbottom sector has to be taken into account, rendering the analysis somewhat more 
involved. We will focus on the scenario with a low value of tan /3 in the following. 

In Fig. ^ we display App^^"^ as a function of the mixing angle for three values of the 
common squark mass, rUq = 150, 250 and 500 GeV, and for tan/3 = 1.6. The contribution is 
practically fiat except for values of the mixing angle very close to the lower limit, 6q ~ — 7r/4. 
This justifies the choice of concentrating on the two extreme cases. In fact, in the case 
Mf^ = Mj:^, the maximal mixing scenario 6q = —tc/A is only obtained exactly when the 
D-terms are set to zero, which is the case when tan/? = 1. One can also have maximal 



^The corrections to the relations between the squark masses will be discussed in subsection 4.4. 

■^As will be discussed below, the case of exact maximal mixing, 9j: = — 7r/4, is not possible in the scenario 
with = and tan P =1.6. Since 6j: is already very close to the maximal value for M/"^ = 200 GeV, 
we will in the following refer to this scenario as "maximal mixing" . 

'^The correction to Ap discussed here directly corresponds to a correction to the effective parameter T 
defined in Ref. |Q, or equivalently to ei as given in Ref. or the combination of parameters defined in 
Ref. j2^. Using the 1997 precision data, the resolution on ei is estimated to be 1.3 x 10^"^ p5|. A similar 
estimate can be readily obtained using the present experimental errors on the world average pq | of Mw , 80 
MeV, and sin^ 6i^, 2.2 x 10"^, in eq. 
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mixing if the sums of and with their corresponding D-terms are equal, making the 
diagonal entries in the mass matrices eqs. (1-2) identical. 

The parameter M^^ docs not only influence the mixing but also the mass splitting 
between the scalar top quarks. The effect of varying M^"^ on Apg^^^ is displayed in Table 1 
for the case of exact maximal mixing, i.e. tan/5 = 1. The scenario with tan/3 = 1.6 yields 
similar numerical results. For large values of M^^ the contribution to the p parameter 
can become huge, exceeding by far the level of experimental observability. The increase of 
Apo^^^ with larger values of Af/"^ is due to the increased mass splitting between the two 
stop masses and m^^^ that is induced for nitM^^ ^ m?. For values of comparable to the 
top quark mass, rriq < 200 GeV, the mass splitting is already large even for small M^^ due 
to the additional term in the stop mass matrix; this is similar to the no-mixing case. 

If the GUT relation M^^ ~ M^^ is relaxed and large values of the mixing parameter 
M/'^ are assumed, the sphtting between the stop and sbottom masses is so large that the 
contribution to the p parameter can become even bigger than in the previously discussed 
scenario (see subsection 4.4 for a more detailed discussion). 



rriq (GeV) 


M^^ (GeV) 


ApSUSY X 10-3 


200 


100 
200 
300 

400 


1.68 
1.35 

0.89 
1.28 


500 


200 
500 
1400 

1550 


0.34 
0.23 
2.04 

4.97 


800 


500 
1500 
2000 
3700 


0.12 
0.07 
0.30 
15.8 



Table 1: Apo^^^^ in units of 10 ^ for several values of niq and Ai^^ (in GeV) and tan/? = 1. 
The values of Ml^^ are chosen such that the corresponding squark masses lie in the experimentally 
allowed range. 
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3. Two— loop calculation 
3.1. Renormalization 

The QCD corrections to the squark contributions to the vector boson self-energies, Fig. ^, 
can be divided into three different classes: the pure scalar diagrams (Fig. ^-c), the gluon ex- 
change diagrams (Fig. ^-j), and the gluino-exchange diagrams (Fig. Qc-n). These diagrams 
have to be supplemented by counterterms for the squark and quark mass renormalization 
(Fig. ^-c) as well as for the renormalization of the squark mixing angle, (Fig. ^). The 
three different sets of contributions together with the respective counterterms are separately 
gauge-invariant and ultraviolet finite. For the gluon exchange contribution we have only 
considered the squark loops, since the gluon exchange in quark loops is just the SM contri- 
bution, yielding the result Apf^ = -Ap™|i^(i + ^rVS) Q. 

As mentioned above, the results presented in the following are precisely the same in 
dimensional regularization as in dimensional reduction. The renormalization procedure is 
performed as follows. We work in the on-shell scheme where the quark and squark masses 
are defined as the real part of the pole of the corresponding propagators. One further needs 
a prescription for the renormalization of the squark mixing angle. The renormalized mixing 
angle can be defined by requiring that the renormalized squark mixing self-energy ir^^^^^{q'^) 
vanishes at a given momentum-transfer q^, for example when one of the two squarks is on- 
shell. This means that the two squark mass eigenstates qi and q2 do not mix but propagate 
independently for this value of g^. Expressing the parameters in eq. (^ by renormalized 
quantities and choosing the field renormalization of the squarks appropriately, this renor- 
malization condition yields for the mixing angle counterterm 

SWo) = 2 ^ 2 %.-.(9g). (14) 
ml - ml 

Finally, we have also included, as a check, the field renormalization constants of the quarks 
and the squarks in our calculation; they of course have to drop out in the final result, which 
we have verified by explicit calculation. 

The one-loop diagrams of Fig. ^-c provide the renormalization of the squark masses 
and the squark wave functions. As the one-loop squark contributions to the vector boson 
self-energies are finite at vanishing external momentum (see eq. (^)), we notice that the 0(e) 
part of the squark mass counterterm is not needed. We have contributions from the three 
diagrams of Fig. b, and c involving gluon exchange, gluino exchange and a pure scalar 
contribution, respectively. The explicit form of the pure scalar contribution is not needed as 
we will see later. The contribution of the gluon exchange to the squark mass counterterm is 
given by 

ma o^mn. = mz 

while the one of the gluino exchange reads 

{ml + rra? - m? ) BQ{m\,mq,mg) + AQ^m-g) + Ao(m,) 
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(15) 



maOmg. = 
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+ 2(-l)* sin26'^ niq BQ{mi^,mq,mg) 



(16) 



where Aq and Bq have been defined previously, and are needed only up to 0(1) in the e 
expansion. Although in our renormalization scheme the contribution of the quartic-squark 
interaction to the squark masses drops out in the final result, we will also give its expression 
for later convenience ii' = ?> — i): 



OTT 



cos^ 2egAoimqJ + sin^ 2egAoimg^, )] . (17) 



Concerning the quark mass counterterm. Fig. |d, in principle the 0(e) term is needed because 
the quark loop contributions to the vector boson self-energies are ultraviolet divergent even 
at = 0. However, in Ap this contribution drops out as the one-loop quark contribution 
to this physical quantity is finite. As mentioned above, the gluon contribution to the quark 
mass counterterm is only relevant for the pure SM correction and is not needed in the present 
context. The gluino contribution can be expressed as 

S'^rug = V (-1)' sin26'g mg Eo(m^,mg,mg,) 

Finally, the counterterm for the squark mixing angle, defined at a given q^, is given by 

^HQo) = \4:mtmgBo{qQ,mt,mg) + sin2^i(Ao(mi ) - Ao(mj ))j . (19) 

As discussed above, for the value of q^ one can either choose m^^ or m^^, the difference being 
very small. In our analysis we have chosen = m?^. This renormalization condition is 
equivalent to the one used in Refs. for scalar quark decays. 



Let us now discuss the separate contributions of the various diagrams. The contribution 
of the pure scalar diagrams vanishes, while for the gluon-exchange diagrams one needs to 
calculate only the first four genuine two-loop diagrams. Fig. ^-g, and the corresponding 
counterterm with the mass renormalization insertion, Fig.|^. The other diagrams do not 
contribute for the following reasons: 

(a) The diagram Fig. ^ is exactly canceled by the corresponding diagram with the mass 
counterterm insertion. Fig. It should be noted that the expression for Bq is needed 
in order to obtain this result. 

(b) The reducible diagram of Fig. ^a involving the quartic squark interaction contributes 
only to the longitudinal components of the vector boson self-energies and can therefore 
be discarded. 
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(c) The diagrams Fig. ^b-c are canceled by the corresponding diagrams with the countert- 
erms, Fig. ^-b, for mass renormahzation (for the diagonal terms) and mixing angle 
renormalization (for the non-diagonal terms). The diagrams in Fig. |^d-f for the vertex 
corrections contain only field renormalization constants which drop out in the final 
result. 

(d) The gluon tadpole-like diagrams of Fig. ^-j give a vanishing contribution in both 
dimensional regularization and reduction. 

For the gluino-exchange diagrams, one has to calculate all diagrams of the types shown 
in Fig. ^-n and their corresponding counterterm diagrams depicted in Fig. ^. 

We now briefly describe the evaluation of the two-loop diagrams. As explained above, we 
have both irreducible two-loop diagrams at zero momentum-transfer and counterterm dia- 
grams. After reducing their tensor structure, they can be decomposed into two-loop scalar 
integrals at zero momentum-transfer (vacuum integrals) and products of one-loop integrals. 
The vacuum integrals are known for arbitrary internal masses and admit a compact repre- 
sentation for e ^ in terms of logarithms and dilogarithms (see for instance Ref. [|2^), while 



the one-loop integrals Aq,Bq are well-known (see eq. (0)). We have used two independent 
implementations of the various steps of this procedure and obtained identical results. 

In the first implementation, the diagrams were generated with the Mathematica package 
FeynArts |^ . The model file contains, besides the SM propagators and vertices, the relevant 



part of the MSSM Lagrangian, i.e. all SUSY propagators (ti, ^2, &2, ^) needed for the 
QCD-corrections and the appropriate vertices (gauge boson-squark vertices, squark-gluon 
and squark-gluino vertices). The program inserts propagators and vertices into the graphs in 
all possible ways and creates the amplitudes including all symmetry factors. The evaluation 
of the two-loop diagrams and counterterms was performed with the Mathematica package 
TwoCalc By means of two-loop tensor integral decompositions it reduces the amplitudes 
to a minimal set of standard scalar integrals, consisting in this case of the basic one-loop 
functions Aq,Bq (the Bq functions originate from the counterterm contributions only) and 
the genuine two-loop function T134 [^, i.e. the two-loop vacuum integral. As a check of 
our calculation, the transversality of the two-loop photon and mixing self-energies at 
arbitrary momentum transfer and the vanishing of their transverse parts at = was 
explicitly verified with TwoCalc. Inserting the explicit expressions for Aq, Bo,Ti34 in the 
result for Ap, the cancellation of the and 1/e poles was checked algebraically, and a 
result in terms of logarithms and dilogarithms was derived. From this output a Fortran code 
was created which allows a fast calculation for a given set of parameters. 

In the second implementation, completely independent, the diagrams were not generated 
automatically, but the analytic simplifications and the expansions in the limiting cases (small 
and large gluino mass, maximal and minimal mixing) were carried out by using the Math- 
ematica package ProcessDiagram In this way the results can be cast into a relatively 



compact form, shown in the following. 
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4. Two— loop results 
4.1. Gluon exchange 

In order to discuss our results, let us first concentrate on the contribution of the gluonic 
corrections and the corresponding counterterms. At the two-loop level, the results for the 
electroweak gauge-boson self-energies at zero momentum-transfer have very simple analyt- 
ical expressions. In the case of an isodoublet {u, d) where general mixing is allowed, the 
structure is similar to eq. @j and eq. (j^) with the gyMj ^ given previously: 



n 
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(0) 



nzz(o) 
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(20) 



q — u,d 
j,j = l,2 



The two-loop function Fi{x,y) is given in terms of dilogarithms by 



xy X 

Fi(x,y) = X + y — 2 m — 

x-y y 



2 + - In - 

y y 



+ 



[X 



yW 



\x 



yy 



X 



y 



2(x-y)Li2 1 



• (21) 



This function is symmetric in the interchange of x and y. As in the case of the one-loop 
function Fq, it vanishes for degenerate masses, Fi{x^x) = 0, while in the case of large mass 
splitting it increases with the heavy scalar quark mass squared: Fi{x, 0) = a;(l + 7r^/3). 

From the previous expressions, the contribution of the (t, b) doublet to the p parameter, 
including the two-loop gluon exchange and pure scalar quark diagrams are obtained straight- 
forwardly. In the case where the b mixing is neglected, the SUSY two-loop contribution is 
given by an expression similar to eq. (|l^): 



SUSY 
, gluon 



+ cos^ e^Fi 



-sin^^t-cos^ OtFi (' 



rrir , m: 



't2 



2 2 

rrir , mr 



(22) 



The two-loop gluonic SUSY contribution to Ap is shown in Fig. |^ as a function of the 
common scalar mass rriq for the two scenarios discussed previously: % = and 9}: ~ — 7r/4. 
As can be seen, the two-loop contribution is of the order of 10 to 15% of the one-loop 
result. Contrary to the SM case (and to many QCD corrections to electroweak processes 



in the SM, see Ref. |^ for a review) where the two-loop correction screens the one-loop 



contribution, Ap^gf^^j^ has the same sign as Apg^^^. For instance, in the case of degenerate 
scalar top quarks with masses ^ m^, the result is the same as the QCD correction to 
the (t, b) contribution in the SM, but with opposite sign, see section 4.3 below. The gluonic 
correction to the contribution of scalar quarks to the p parameter will therefore enhance the 
sensitivity in the search of the virtual effects of scalar quarks in high-precision electroweak 
measurements. The dependence of the two-loop gluonic contribution on the stop mixing 
angle % exhibits the same behavior as the one-loop correction: Ap^^p^^^ is nearly constant 
for all possible values of 6f, only in the region of maximal mixing it decreases rapidly. 
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4.2. Gluino exchange 

Like for the gluon-exchange contribution we have also derived a complete analytic result 
for the gluino-exchange contribution. The complete result is however very lengthy and we 
therefore present here explicit expressions only for the limiting cases of light and heavy gluino 
mass, and for the case of no squark mixing. The complete expression is available in Fortran 

and Mathcmatica format from the authors. 

In order to make our expressions as compact as possible, we use = sin 6f and = cos 6^ 
as abbreviations and introduce the following notation, 

dxy = {ml - ml) , where mi = m^^, m2 = m^^, = m^^, = m^. 

The gluino-contribution for vanishing gluino mass is given by 

as CfNcGf 



SUSY 
uino 



TT 




16v/27r2 




o 2 2 J // 2 I 2\r 2,2 2 2\ 4|22\A 

- 2sj-m^^rfi2((m^^ + mt){mi^ +mt- m^^ - m^J - + m^mfj j 

- o?i2(2m|(m?^ + m^^ - m|) 

+ 2m?^m^(m?^m?^ - + mf) - m^im^ + 3mf^)) + s}(mj^{ml^ - mf^d 
+ 2m^m| (m|mf^ (mf^ - m| + 3m|) - 2m|(m| + J + mfj 
+ mtiSml^ + ml)d2i + 2m^(4i + mf^dis)) 

+ '^4^l{dli{m\m\ - ml 

I 2/4/4 4 \j io22/2 2/2 2 2\i 6\ 

+ - ^h)d2i + 2m,m,-^(m,-^m^Jm^^ - m^^ - mQ + m^-J 

+ m^(3m|^ + m|)(ii2 + Im^i-dl^ + mf^ci2i)) 

- 24^? (42(^1^1 - ^t))) 

In^fm- ) / 
+ ,0 'y {dU2mt mj + 2m? 

- ^ti^^^L + + - '^s\ml^di2{{ml^{ml^ + m| + m|) 

o2 2 2\, 2/2/2, 2 o2\, 2 2\, 4/o 2 2 2 \\\ 

- ^mim^^mQ + {mi^{mi^ + m^-^ - 3mgJ + m^-^mj-J + (2m^^ - m^-^ - mQ) j 
-21n(mfjln(m|) ^ ^ 2 , -21n(mfjln(m|) ^ ^ ^ 

+ ^b^^tiCf + ^ft^^f^^t 
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- 2s}ml^di2{ml{ml^ - m\ - m|) + m? m|)) 

2M<)MMul) 4 .2 , 21n(mgJln(M,,|) ^ ^ 

777-- lllfTTT'- 1 



c| In^ (m| ) 



+ ^-^.^ — T-^i — ;^ -{m%di2{mlm\ + - 2ml m]) 



m}mldLidi2 
+ 2sfdLi(m| + m|)(2m|m| - m? m^^ - rnlrn^^)) 
, 2<ln(m|)ln(Mi,|)cf ^ ^ _ ^ 

^ 2di,ln(m|)ln(|d2t|) ^ 2 4 
ln(m^ ) / 4 

o2 2/ 4/2 2 2\| 2 22, 2rr,r 2 , 2 \ j o 4 

- 2sfmt{m^^{m~^^ - m^^ - m^J + m-^^m^m^^ + mt{2{mf^ + m^JdLi - 2mf^ 
+ 2^lO) + HrnldLi{ml (m? + m|) + 2m?(2m? + m|))) 

+ ^^|^(2m|di^4L + mld,2{^ml + m|) 

- 2m^(42 + ml^dn) - 28'ldl^{m}m}^ - m^)) 
mg5fln(m|)ln(mg) ^ 2 . 2 2 . 2 ^ 2 x ^ 2 /o 2 02 .2 

i 1 ^2 

+ 2m^((m| + mlf - m^iSml + m|))) - 2sfdL2(m 
{2ml m}^ - ^ti'ml + ml))) 



^ 2m|<4c fln(m|)ln(|du|) 



24s|ln(mnin(|rf2*|) 



mf diiiml di2 + mf (iia) - "^i 



14 



{dL2{m\ + m|) + m^(6m?^ - 2m| - 4m|)) 
- 2s^^dL2ml{m\{m\ + m|) + 2m,2(m| + 2m| 
In(m^) 



mj mj dLidL2\^ ^ 



+ m^{6ml^{ml^ - m| - m|) + 2m|(m| + 2m|))) 

- sf[dLidL2m^^{10mlml - 3m,2(m| + m|))^ 

- 4s?m2(-c/2^ + 42 - - 2m| + m|)) + 2sfm'idLi{ml - 3m|; 

I '^2tl^(M2t|) / 2 J J A 2 2 2 J io4 2j/o2 2- 

+ 9^2 ^2^ ^ V^hdL2di2 - As^m^m- dLi + 2s^m^dL2{3m- - 

^^^2 t L2"12 ^ 



- 6m^) + s?(6m2(-3m| + m| - Gm^) + 7r2((m| + m?^)^ - (m| + m?^)^ 

- 2m2rfi2)) + Qsfmliml + m| + 6m?)) . (23) 

The result for rrig = is compared to the complete result of the gluino contribution for 
rrig = 10, 200, 500 GeV in Figs. ||, In Fig. ^ the results for the no mixing scenario are 
displayed, while Fig. |^ shows the maximal mixing case. As expected, the curves for rUg = 
and rrig = 10 GeV are very close, while the curves for large gluino masses significantly 
deviate from the rrig = case. For the no mixing scenario the light gluino expression of 
eq. (p3[) reproduces the exact result within 5xl0~^ for rrig < 90, 130, 350 GeV in the cases of 
rrig = 100, 250, 500 GeV, respectively. In the maximal mixing scenario the correction is much 
smaller for a light gluino than in the no mixing case; the light gluino expression reproduces 
the exact result within 5x 10"^ for m-g < 20, 300, 400 GeV in the cases of rrig = 100, 250, 500 
GeV and maximal mixing. 

As second limiting case we give a series expansion in powers of the inverse gluino mass 
up to C(l/m?). Note that, due to the term linear in rrig in eq. ([T8|), the expansion actually 
starts at 0{l/mg). 
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Fig. shows the quahty of the heavy gluino expansion up to 0(l/m?) in the two cases 
of minimal and maximal mixing for m^ = 100 GeV. The expansion of eq. P^ ) reproduces 
the exact results within at most 30% (which corresponds to a maximum deviation of about 
5 X 10~^ in Ap) for m^ > 200 GeV in the no mixing case and for m^ > 340 GeV in 
the maximal mixing scenario. We have also verified that higher orders in the expansion 
improve significantly the convergence of the series, but we have not included unnecessary 
long expressions. As can be expected, when the value of m^ is increased, the quality of the 
heavy gluino approximation deteriorates. At the same time, however, the gluino correction 
becomes smaller and can mostly be neglected; it never exceeds 1.5 xlO~'^ and lxlO~"^ for 
mq > 200 GeV in the cases of minimal and maximal mixing, respectively (see Figs. ^ and 
1^). We conclude that in the region of the parameter space where the gluino contribution 
is relevant, and the gluino is not too light (say m^ > 300), eq. ( plj) approximates the full 
expression sufficiently well. 

As mentioned above, for most of the parameter space the mixing in the stop sector is 
either zero or nearly maximal. As a third limiting case we give Ap^gf^fj^^ for arbitrary m^, 
but with zero stop mixing. The result in this limiting case has already been displayed in 
Fig. H and Fig. (solid curve). 
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Here we have used X{x,y, z) = y{z — x — y)'^ — Axy and the function $(x, 2;) is defined 
according to the sign of X{x, y, z)"^. For A(x, z)"^ > 0, 
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(26) 



(27) 



where Cl2(a;) = ImLi2(e*^) is the Clausen function; limiting expressions and additional 



details on this function can be found in Ref. 28 . 



4.3. Supersymmetric limit 

The sign of the gluonic two-loop contributions is, as discussed above, always the same as the 
sign of the one-loop contribution, contrary to the case of the two-loop contributions in the 
SM. In the limit of vanishing gluino mass, mi^ = m^^ = mt, m^^ = = (supersymmetric 
limit), the gluon exchange contribution of the scalar quarks reads 

ApSn = "^^^"^^'^'^ = Apl^^-^ (1 + -], (28) 
Pi.giuon ^27^3 12 3 vr V 3 y ' ^ ' 
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which exactly cancels the quark loop contribution |T^. The gluino exchange contribution in 
this limit is given by 



^Pl,gluino - -^Po 3—' l^yj 

which numerically cancels almost completely the contribution of eq. (^). 
4.4. Corrections to the squark masses 

The analytical formulas for Ap given in the previous sections are exclusively expressed in 
terms of the physical squark masses. The formulas therefore allow a general analysis since 
they do not rely on any specific model assumption for the mass values. For our numerical 
analysis, however, since the physical values of the squark masses are unknown, we had to 
calculate the masses from the (unphysical) soft-SUSY breaking mass parameters. We have 
concentrated on the MSSM scenario with soft-breaking terms obeying the SU(2) relation 
Mi^ = . The t and h masses are obtained by diagonalizing the tree-level mass matrices. 
Using eqs. (|T]-§), the masses m^^, m^^ m^^ = m^^Q as well as the mixing angle 9i can be 
expressed in terms of the three parameters M^^, Mj^, and M^^. This yields the tree-level 
mass relation 



= cos^ 6'tm| + sin^ 6'jm| - mf - cos 2/3 M^. (30) 

When higher-order corrections are included, the quantities M^^ and M^^ are renormal- 
ized by different counterterms once the on-shell renormalization for the t and b squarks is 
performed. Requiring the SU(2) relation 

M? +5M? =Mi +5 Mi (31) 

for the bare parameters at the one-loop level, one obtains Mi^ 7^ M-^^ for the renormalized 
parameters (see also the discussion in Ref. ^^). The difference ^M"^^ — (5M| = Am? with 

5Mj^ = cos^ 9iSm\ + sin^ 6'i5m| + (m| - m| ) sin 29^59^ - 2mt5mt (32) 
5Ml = 5ml (33) 



constitutes a finite 0{as) contribution to the 61 mass compared to its tree-level value. The 



mass shift Am? can also be obtained by replacing the tree-level quantities in eq. (pOD by their 
renormalized values and the corresponding counterterms. We have explicitly checked that 
the quantity Am? is indeed finite. Note that the squark mass and mixing angle counterterms 
in eq. (p^-pBf) also receive contributions from the pure scalar diagrams of the type of Fig. ^ 
which in the two-loop results given above canceled between the counterterm graphs and 
the two-loop diagrams. For the top mass counterterm in eq. ( p3D also the graph with 
gluon exchange enters, which was absent in the results given above. Due to the latter 
contribution the 0{as) correction to eq. ( pUD is different in the dimensional regularization 

^Since we set mi, = 0, thus neglecting mixing in the b sector, we have = m-^i,' ^^"^ "^bj decouples 
from the p parameter. 
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and the dimensional reduction scheme, while for all results given above the two schemes 
yield identical results. In our numerical analysis in this section we use the result for the top 
mass counterterm in dimensional reduction. 

For the two-loop contribution to the p-parameter, the difference generated be using 
the tree-level mass relations instead of the one-loop corrected ones is of three-loop order 
and therefore not relevant at the order of perturbation under consideration. However, if 
the one-loop contribution to Ap is evaluated by calculating the squark masses from the soft 
breaking parameters as described above, the 0{as) correction to the bi mass should be taken 
into account. This gives rise to an extra contribution compared to the results discussed in 
section 2. 

The tree-level relation M^^ = M^^, which has been widely used in our numerical exam- 
ples, can in principle be maintained also for the renormalized parameters at the one-loop 
level with on-shell ti,t2 mass renormalization. Alternatively, in the spirit of the discussion 
given above, one may assume that the symmetry is extended to the bare parameters: 

M? + 6Mf =Mj + 5Mj . (34) 

th th iR ^ ' 

Accordingly, the squark masses and the mixing angle are in this scenario given in terms of 
the two parameters irtq and M^^, and there exist two relations between the squark masses 
and the stop mixing angle at tree-level. 



C2 



L_^cos2,3M| (cf, - =in^9, (| - . (36) 



*i *i * sm^ 61,- - cos2 6^- 



At 0{as), the right-hand sides of eqs. (|35|) and (^) receive the extra contributions Am? 



t2 

and Am^ , respectively, where 



, ICSpCULlVCiJ', WiiClC 

sin 29r6er 



Ami = - - 2mt6mt + — — f'^'^^t^^t cog 2/? Ml(-1 + -sL). (38) 

''1 *i ^1 2(sin2ei-cos2^^)2 ^ 3 ly; V ; 



These relations have been used in the numerical evaluation displayed in Fig. |TT]. The differ- 
ence to the case with M^^ = M^^ at the renormalized level is only marginal. 



The effect of the 0{as) corrections to the mass relations on Apo is shown in Fig. \L1 
where the one-loop correction to the p parameter is expressed in terms of rriq, but with the 
corrections eqs. (|37|-p^) to the mass relations taken into account. The results are shown for 
the no-mixing and the maximal-mixing case and for two values of the gluino mass, rrig = 200 
and 500 GeV. Compared to Fig^, where the tree-level mass relations have been used as 
input, the contribution of the {t, b) doublet to the p parameter is reduced for the parameter 
space chosen in the figure. 

The inclusion of the one-loop corrections to the squark mass relations does not always 
lead to a decrease of Apo, but can also give rise to a significant enhancement. This is 
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quantitatively shown in Table 2 for the same scenario as in Table 1, i.e. for the maximal- 
mixing case with tan/5 = 1. The shift for the sbottom mass in this case follows from 
eqs. (|32|-p3D in the limit 9i = —it /A with 66j: = 0. The numerical results are very similar to 
the case with tan/? = 1.6. The value chosen for the gluino mass is rrig = 500 GeV, and for 
completeness also the full two-loop contribution Apf^^^ is given. One can see that for large 
values of the non-diagonal element in the t mass matrix, Apo can become larger compared 
to the entries in Table 1, thus significantly increasing Ap at the two-loop level in a range 
where the one-loop contribution is already quite large. 



rrig (GeV) 


M,^^ (GeV) 


ApSUSY X 10-3 


ApSon X 10-^ 


A SUSY X in-4 
'-^Pl.gluino ^ -"-U 


200 


100 


1.17 


1.92 


0.09 




200 


0.90 


1.53 


0.61 




300 


0.57 


1.05 


1.37 




400 


1.22 


1.71 


2.80 


500 


200 


0.23 


0.39 


-0.08 




500 


0.13 


0.27 


0.52 




1400 


2.50 


2.68 


3.94 




1550 


5.71 


6.18 


5.34 


800 


500 


0.07 


0.14 


0.13 




1500 


0.07 


0.11 


0.99 




2000 


0.31 


0.42 


2.31 




3700 


14.23 


16.89 


16.05 



Table 2: Apg^SY ^^^^^ ^q-s f^^. ^.j^g 

same scenario as in Table 1, i.e. with tan/3 = 1, for 
several values of niq and M^^. The 0(a<j) correction to the tree-level mass relation has been taken 
into account for rrig = 500 GeV. The two-loop contributions Ap^ gp^^^^ and Ap^ gf^^^^ are also given 
(in units of 10-^, rug = 500 GeV). 



In the situation where the relation M^^ = M^^ is relaxed and the squark masses and the 
mixing angle are derived from Mj^, M^^ and M^^, Fig. [l^ shows Apg^^^ for the two choices 
MijMi^ = 300/1000 and 1000/300 as a function of M/^^. For the MijMi^ = 300/1000 
case the solid line corresponds to the use of the tree-level masses, while the others reflect 
the use of the one-loop correction Am? to the tree-level masses following from eq. (|3^^3|) 
for two gluino masses, nig = 200,500 GeV. For the M^^/M^^ = 1000/300 case the result is 
given by the dotted line. It is insensitive to the one-loop correction to the squark masses; 
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we therefore show only a single curve. For the M^^/Mf^ = 300/1000 case, Apo is decreased 
for small values of M^^ but is increased for large M^^. The effect is more pronounced for 
heavier gluinos. For large values of M^^^ Apo can become huge in this scenario, exceeding 
the level of experimental observability. The bounds on the breaking parameters obtained 
from experiment will therefore crucially depend on the proper inclusion of the two-loop 
contributions. 

5. Conclusions 

We have calculated the 0{as) correction to the squark loop contributions to the p parameter 
in the MSSM. The result can be divided into a gluonic contribution, which is typically of 
(9(10%) and dominates in most of the parameter space, and a gluino contribution, which 
goes to zero for large gluino masses as a consequence of decoupling. Only for gluino and 
stop/sbottom masses close to their lower experimental bounds, the gluino contribution be- 
comes comparable to the gluon correction. In this case, the gluon and gluino contributions 
add up to ~ 30% of the one-loop value for maximal mixing. In general the sum of gluonic 
and gluino corrections enters with the same sign as the one-loop contribution. It thus leads 
to an enhancement of the one-loop contribution (expressed in terms of the physical squark 
masses) and an increased sensitivity in the search for scalar quarks through their virtual 
effects in high-precision electroweak observables. This is in contrast to what happens in the 
SM, where the two-loop QCD corrections enter with opposite sign and screen the one-loop 
result. 

While the gluonic contribution can be presented in a very compact form, the complete 
analytical result for the gluino correction is very lengthy. We have therefore not written 
it out explicitly but have given expressions for three limiting cases, namely the result for 
zero squark mixing, for vanishing gluino mass, and an expansion for a heavy gluino mass. 
These limiting cases approximate the exact result sufficiently well for practical purposes. The 
results have been given in terms of the on-shell squark masses and are therefore independent 
of any specific scenario assumed for the mass values. For different scenarios we have analyzed 
the extra contributions caused by the 0{as) correction to the tree-level mass relations. 
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5i cb 

V /' ~\ V 

vAAAAAaW wwwv^ 

V \^ /' V 

Figure 1: Feynman diagrams for the contribution of scalar quark loops to the gauge boson 
self-energies at one-loop order. 
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Figure 2: One-loop contribution of the (t, b) doublet to Ap as a function of the 
common squark mass rriq for Of = and 6f n /4 (with tan/? = 1.6, M^^^ = 0, and 
M/'^ = or 200 GeV). 
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Figure 3: Dependence of the one-loop contribution Apg'^^'^ on the stop mixing 
angle 9^. The parameters are the same as in Fig. |^. 
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V ,'''*"'\ /''*^''\ V V / 4 9k\ V V \ / V 

vAyWVW • WWVW vAAAAAA/i tAAAAA/V> vAAAAAAAAATVW^VWVWW^ 

(a) (b) (c) 



vTVAAAAAA ^ WVWNA/- vAAAAAAATRTOOTWipkAAAAA/V' vAAAAAAA WWVW 

C^' '57' 'X 

(d) (e) (f) 



V / '\ V V \^ V V / 3i qi'', V 

^A/WWvriyiS^ WWWV^ ^/VWWW^A7^y^AAAAAAAAA/^ ^AAAA/\Ay> WV\A/\A/^ 




Figure 4: Feynman diagrams for the contribution of scalar quark loops to the gauge-boson 
self-energies at two-loop order. 



29 




(a) (b) (c) 




(d) (e) (f) 

Figure 5: Counterterm contributions to the gauge-boson self-energies at two-loop order. 
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Figure 6: One-loop diagrams contributing to the squark mass and mixing angle countert- 
erms. 
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Figure 8: Ap^gf^fj^^ as a function of rriq in the no-mixing scenario; tan/3 = 1.6 and 
nig = 0, 10 (the plots are indistinguishable), 200 and 500 GeV. 
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Fi gure 9: ^Pi^uifno function of Triq in the maximal mixing scenario; tan j3 
and m-g = 0, 10, 200, and 500 GeV. 
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Figure 10: Comparison between the exact result for Ap^gjjj^^ and the result of the 
expansion eq. ( pi]) up to 0(l/m|) for the two scenarios of Fig. Q niq = 100 GeV. 
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Figure 11: Apg^^^ calculated with one-loop corrected squark masses as a function of 
rriq for Oi = and 9i ~ -7r/4 (tan (5 = 1.6, M/"^ = or 200 GeV) and two values of m~g. 
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Figure 12: Apg^SY ^ ^ function of M^^ for tan/3 = 1.6 and for MijMi^ = 1000/300 
(dotted line, the lines for tree-level and one-loop parameters are not distinguishable) and 
300/1000 (solid line for tree-level parameters, for the one-loop parameters: dash-dotted 
for rUg = 200 GeV and dashed for rrig = 500 GeV. 
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